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Estimates of the quantum accuracy threshold often tacitly assume that it is possible to interact 
arbitrary pairs of qubits in a quantum computer with a failure rate that is independent of the distance 
between them. None of the many physical systems that are candidates for quantum computing 
possess this property. Here we study the performance of a concatenated error-detection code in a 
system that permits only nearest-neighbor interactions in one dimension. We make use of a new 
message-passing scheme that maximizes the number of errors that can be reliably corrected by the 
code. Our numerical results indicate that arbitrarily accurate universal quantum computation is 
possible if the probability of failure of each elementary physical operation is below approximately 
10~^. This threshold is three orders of magnitude lower than the highest known. 
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I. INTRODUCTION 

For a quantum computer to reliably outperform a clas- 
sical computer, it must be robust against the effects of 
decoherence and imprecise quantum control. There are 
many ideas on how to achieve such fault tolerance, in- 
cludingtopological quantum computing 0, H, 0] , surface 
codes III, 0, m , color codes (H , self-correcting codes , 
and concatenated codes [1, The threshold theorem in- 
dicates that, under certain conditions, scalable quantum 
computing is possible in principle [TO, [Hi • The theorem 
asserts that if the probability of failure of each elemen- 
tary physical operation in a quantum computer is below 
some threshold then arbitrarily accurate quantum com- 
putation can be performed efficiently given sufficient time 
and qubits. The actual value of the threshold for a given 
error-correction code depends on a number of parameters 
that describe the quantum computer and the noise that 
affects it [12]. 

In most estimates of the accuracy threshold, the as- 
sumption is made that it is possible to interact arbitrary 
pairs of qubits in a quantum computer with a failure rate 
that is independent of the distance between them. This 
property is desirable, since it allows higher failure rates 
to be tolerated, but it is unrealistic in the limit of many 
qubits. In all of the physical systems that are candidates 
for quantum computing the range of controllable inter- 
actions is constrained such that at least some pairs of 
qubits will need to be transported before they can un- 
dergo logic gates such as cnot and cphase. Moreover, 
we expect that in many systems, especially trapped-ion, 
superconducting, and solid-state systems, hardware lim- 
itations will require that the qubit array is two- or one- 
dimensional [il [ii, [H, [3 . 

The combination of local interactions and low coor- 
dination number is not a problem in principle as it is 



known that the threshold theorem still applies [TO, [13] • 
The value of the threshold has been quantified for sys- 
tems that permit only nearest-neighbor interactions in 
two dimensions [H, [TqI, [2Q|| and in various quasi one- 
dimensional settings [2l|, [2211, a nd the threshold is also 
known for superconducting and ion-trap [1^ archi- 
tectures with similar properties. However, there are a 
large number of systems under development that per- 
mit only nearest-neighbor interactions in one dimension 

[n, m, m, 113, m, h iio, [m, [32, s, s, [si, iip, for 

which the threshold is unknown. The threshold in one 
dimension is expected to be significantly lower than in 
all other cases. 

Here we find the accuracy threshold for a system that 
permits only nearest-neighbor interactions in one dimen- 
sion. Where required, qubits are transported via nearest- 
neighbor SWAP gates. To minimize this overhead we use 
a small error-detection code. To correct errors we make 
use of a new message-passing method that uses classi- 
cal information gathered during error detection to maxi- 
mize the number of errors that can be reliably corrected 
(37! . Our numerical results indicate that arbitrarily accu- 
rate quantum computation is possible if the probability 
of failure of each elementary physical operation is below 
approximately 10~^. 



II. [[4,1,2]] SUBSYSTEM CODE 

The code that we have chosen to use is the [[4,1,2]] 
subsystem code, a stabilizer CSS quantum code that en- 
codes one logical qubit into four physical qubits [3, [l^ • 
Its stabilizer, 5, is 
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where Xi and Zi represent the Pauli operators (Jx and (Jz 
applied to the i^^ qubit respectively. Identity operators 
and tensor products between operators are omitted. Al- 
though there are two degrees of freedom in which to store 
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encoded information, we choose to ignore one of these en- 
coded qubits. Then, elements in the non-Abelian group, 
T, generated by the operators 

X1X2 , X3X4 , Zi Z3 , Z2 , (2) 

act trivially on the sole protected qubit. The encoded Z 
and X operators for this qubit are 

XiX3,ZiZ2. (3) 

Because the elements in S are products of elements in T, 
and since all elements in T commute with all elements 
in S and the encoded operators, to determine the eigen- 
values of the elements of S it is sufficient to measure the 
eigenvalues of the elements in T . This may change the 
state of the system but it will not affect the state of the 
protected qubit. This property allows us to use only two 
ancilla qubits to simultaneously perform operator mea- 
surements of the operators in T [39]. We decode the 
syndrome by taking the parity of the two measurement 
outcomes in each basis. This is used to infer the pres- 
ence or absence of errors in each basis — if the parity is 
even there is no error, if the parity is odd there is an 
error. For the [[4,1,2]] subsystem code, cnot is a valid 
transversal gate and is a valid transversal gate up to 
a permutation of qubits. 

III. USING CONCATENATED ERROR 
DETECTION TO CORRECT ERRORS 

In our scheme the [[4,1,2]] subsystem code is concate- 
nated such that physical qubits form encoded qubits 
which in turn form higher-level encoded qubits and so 
on. However, while concatenating the code / times re- 
sults in a final code that has distance 2^+^, if each level 
of error detection operates independently of every other 
level, then the code cannot reliably correct even a single 
error — a single physical failure may cause an encoded er- 
ror at any level. To do more than simply detect errors 
classical messages must be passed from each level to the 
level above. These messages serve to indicate the loca- 
tion of potential errors, thereby removing the ambiguity 
in the cause of any odd parity syndrome that is observed, 
which allows errors that are detected to be located and 
corrected. 

The message-passing method applied in this paper is 
similar to those of Refs. [10, [U, [42[ but, unlike those, it 
will only fail if the number of concurrent errors is greater 
than or equal to half of the distance of the final con- 
catenated code — in the case of the [[4,1,2]] subsystem 
code, the leading-order exponent of the probability of 
failure scales with the number of levels of concatenation 
as l(physical), 1, 2, 4, 8, 16, 32, and so on. This is also 
a property of the noisy-channel method of Ref. [i^ . Al- 
though we will only describe the message-passing method 
in the case where it is applied to the [[4,1,2]] subsys- 
tem code, it can be applied to codes with greater dis- 



tance. A detailed description of the method is contained 
in Ref. 13- 

Messages in the method consist solely of classical in- 
formation which indicates our confidence that each lo- 
cation in the circuit has not failed given what is known 
about it. We call these messages flags. The probabil- 
ity that any given flag represents an actual error is de- 
scribed by the weight, carried by the flag such that 
Pr (error) = O(p^), where p is the probability of failure 
of an elementary physical location in the quantum com- 
puter. Since we assume that errors at each location are 
independent of all other locations, the probability of a 
set of flags representing actual errors is 

C)(pE^), where 

the sum is over the entire set. Flags are raised at every 
elementary physical and encoded location. All flags at 
the physical level are given weight equal to one by defini- 
tion and at all other levels weights are determined during 
error detection at the level below. Flags in the X and Z 
bases are separate. 

As flags are raised they are classically propagated 
through the error-detection circuit to determine the effect 
that an error at the location at which the flag originated 
would have on the data and ancilla qubits at the point 
of syndrome extraction. Note that cnot copies X errors 
from the control qubit to the target qubit and Z errors 
in the opposite direction, and that we define the point of 
syndrome extraction to be immediately after the ancilla 
qubits are measured. Since there are significantly fewer 
effects than there are locations at which flags are raised, 
flags are binned by effect. Each bin is assigned a single 
weight which is only updated if a flag is raised with a 
weight lower than its existing value. 

Figure [1] illustrates the binning system where, for sim- 
plicity, non-local circuits are used and only X errors are 
considered. Note that the X-error syndrome is obtained 
using the second half of the circuit in Fig. [H so we are 
interested in the effect of X errors at the second set of 
measurements. X errors do not affect the Z-error syn- 
drome, but since they can occur during the circuit to ob- 
tain the Z-error syndrome this part of the circuit must 
be considered. Because of the degeneracy of the code, we 
need only consider two pairs of data qubits rather that 
four individual data qubits. For example, since X1X2 
acts trivially on the encoded qubit, both X\ and X2 give 
the same syndrome and can be corrected by the operator 

Since the failure of a single cnot or SWAP gate can in- 
troduce a pair of correlated errors, these gates require 
special treatment. The weights of the two flags that 
are raised immediately following each two-qubit gate are 
used to update the bin which corresponds to the two- 
qubit correlated error. Note that both of the bins corre- 
sponding to single errors following the two-qubit gate are 
also updated. The weight of the correlated error is taken 
to be the maximum of the two single-qubit weights, as 
opposed to the sum, which would describe the probability 
of the pair of errors occurring without correlation. 

Once the syndrome is measured its most likely cause 
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FIG. 1: Flags are binned according to the effect that an error at the location at which the flag originated would have on the 
data and ancilla qubits at the point of syndrome extraction. For simplicity, only X errors are considered. An X error in the 
bin AGl would change the parity of the ancilla while remaining on one of the first pair of data qubits, an X error in AG2 
would change the parity of the ancilla while remaining on one of the second pair of data qubits, and an X error in A would 
change the parity of the ancilla. Errors in Gl and G2 cannot affect parity of the ancilla, but the weights of Gl and G2 are 
recorded and used to update the weights of AGl and AG2 for the next error-correction cycle. Binning of flags corresponding 
to Z errors is undertaken similarly. 



is identified by finding the bin that has the lowest weight 
while still being consistent with, or matching, the syn- 
drome. In the case of the [[4,1,2]] subsystem code, in 
each error-correction cycle we consider only three bins, 
AGl, AG2, and A. The most likely cause of an odd syn- 
drome is a single error occurring in the bin with the low- 
est weight, as each of these bins corresponds to a change 
in the parity of the ancilla. The most likely cause of an 
even syndrome is always no error at all. If the match 
implies an error on a data qubit then the appropriate 
correction is applied. 

Although we correct for the most likely error it is 
possible that the true error is the complement of the 
correction — that is, the true error and the correction that 
we apply combine to form an encoded operator. To de- 
termine the weight of the flag that is raised in the error- 
correction circuit at the level above, what we will call the 
encoded weight, we calculate the difference between the 
weight of the match on which we act and the weight of its 
complement. For example, if we have an odd syndrome 
and the AGl is the match, then AG2 is the complement 
match and the encoded weight is the difference between 
the weights of the these two matches, AG2 — AGl. 

There is the possibility that a correction will result in 
a state that is outside the code space so that it is neither 
correct nor affected by an encoded operator. Just as the 
complement match is used to determine the probability 
of an encoded error, the conditional probability of single 
errors on the pairs of data qubits can be updated by 
considering other matches. At the end of each error- 
correction cycle the weights of AGl and AG2 are updated 
to the minimum of the weight of any previous locations 
that have not yet had an opportunity to affect the ancilla, 
given by Gl and G2, and a weight obtained during the 
preceding circuit, calculated in an analogous way to the 
encoded weight based on the syndrome that is observed. 

Table I shows a list of all possible flag matches along 
with the various flag-weight updates that would result 



from each of them being acted on, where we have retained 
the notation of Fig. [TJ By the careful consideration of 
every possible cause of every syndrome in this way we 
always have accurate weights for every qubit at every 
level. We can always, therefore, apply corrections based 
on the most likely set of errors. 



IV. ACCURACY THRESHOLD IN ONE 
DIMENSION 

To estimate the threshold for universal quantum com- 
putation in one dimension we construct circuits for error 
detection and for the encoded operations cnot, swap, 
Hadamard, state preparation, and measurement. These 
operations, which we will refer to as CSS operations, 
are sufficient to concatenate error detection and to per- 
form state distillation following the ideas presented in 
Refs. 07 113 • State distillation involves preparing the 
ancillary state = |0) + to enable the logi- 

cal phase gate, and preparing the ancillary state 
iToffoli) = |000) + |010) + 1 100) + I 111) to enable the 
logical Toffoli gate. These gates together with the CSS 
operations complete a universal set for quantum compu- 
tation. Accurate states can be distilled from many noisy 
states provided the noisy states can be made with a fail- 
ure rate lower than some distillation threshold, which is 
typically above one percent. We determine the threshold 
for CSS operations by numerically simulating the circuit 
for an error-corrected logical cnot under a stochastic er- 
ror model. The cnot is chosen because it has the high- 
est failure rate of the CSS operations. As the threshold 
we find is well below the distillation threshold it is the 
threshold for universal computation under our scheme. 

Figure [2] shows a circuit for syndrome extraction for 
the [[4,1,2]] subsystem code where only nearest-neighbor 
interactions on a linear array are permitted. Note that 
the encoded Hadamard can be achieved by transversal 
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None 
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AGl 
AG2 
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AGl + AG2 + A 


Xi/Zi 
none 


A- AGl 
AGl + A 
AGl- A 


A-AG2 
AG2-A 
n/a 


AG2 - AGl 
AGl - AG2 
AGl + AG2 



TABLE I: Table of all flag matches for each of the two possible syndromes [37\- The first column lists the syndrome. The 
second column lists the possible flag matches for each syndrome. The correction and flag update rules, shown in the last four 
columns, depend on which is the lowest-weight flag match. Gi and G2 are weights for Gl and G2 at the current level of error 
correction, which become AGl and AG2 respectively for the next error-correction cycle. Ce is the encoded weight for the 
next- highest level of error correction. 
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FIG. 2: Syndrome-extraction circuit for the [[4,1,2]] subsys- 
tem code on a linear nearest-neighbor array. Q denotes the 
configuration of the encoded qubit after syndrome extraction, 
bold lines indicate data qubits, and plain lines indicate ancilla 
qubits. 



application of the Hadamard gate in addition to remov- 
ing a single SWAP gate from this circuit to permute the 
qubits. This syndrome-extraction circuit has the same 
depth as the non-local circuit in Fig. [1] and the two cir- 
cuits differ only by a rearrangement of qubits and the 
addition of two SWAP gates. This difference is signifi- 
cant, however, since each of the additional SWAP gates 
involves two data qubits. New pairs of correlated er- 
rors are introduced with probability where p is the 
probability of failure of an elementary physical location. 
Without the SWAP gates these particular pairs of errors 
occur with probability 0{p^). Since these errors include 
the encoded operators, this means that an undetected 
encoded error occurs with probability 0{p). 

Figure [3] shows the encoded cnot circuit. A pair of 
correlated errors caused by the failure of one of the SWAP 
gates prior to the transversal cnot results in a pair of 
errors on the encoded control qubit and a pair of errors on 
the encoded target qubit. This is because the transversal 
CNOT copies errors. This means that a pair of correlated, 
undetected encoded errors occurs with probability 0{p). 
The encoded SWAP gate, in which the transversal cnot 
is replaced by a transversal SWAP, does not possess this 
property. Two gates must fail during the encoded SWAP 
for both encoded qubits to fail undetected. 
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FIG. 3: Encoded CNOT circuit for the [[4,1,2]] subsystem code 
on a linear nearest-neighbor array, where notation is consis- 
tent with that of Fig. 2. Replacing all cnot gates with SWAP 
gates gives the encoded SWAP circuit. 



Here we make an important observation: to leading 
order, the relative probabilities of undetected failure of 
first-level encoded locations in the one-dimensional case 
mimic the relative probabilities of failure of physical lo- 
cations in the non-local case. Specifically, all single-qubit 
locations fail with probability a pair of correlated 

errors after a cnot occurs with probability 0{p)^ and a 
pair of correlated errors after a SWAP occurs with proba- 
bility O(p^). A pair of correlated errors can occur after a 
SWAP gate because of a single fault, but they will always 
be detected. A corollary of this is that after the first level 
of concatenation our linear scheme will mimic a non-local 
scheme. In the non-local case we expect to succeed when- 
ever the number of concurrent errors is less than half of 
the distance of the final concatenated code. This implies 
that in the one-dimensional case the first-order exponent 
of the probability of failure should scale with the number 
of levels of concatenation as l(physical), 1, 1, 2, 4, 8, 16, 
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and so on. 

Note that Figs. [2] and [3] show circuits constructed from 
physical gates. To generate circuits for error correction 
and encoded operations at higher encoded levels we re- 
place all physical gates with encoded gates. Like in other 
concatenated schemes, after every encoded operation we 
perform error correction on the encoded qubits involved 
in that operation All circuits are designed so that 
ancilla qubits are always available to perform error cor- 
rection using the circuit in Fig. [2l 

To attempt to verify that our scheme performs as ex- 
pected and to determine the threshold, we simulate a 
logical CNOT under a stochastic error model. The simu- 
lated circuit consists of a CNOT extended rectangle pjj . 
the failure rate of which is meant to approximate the 
failure rate of a CNOT in some algorithm or at a higher 
level of concatenation. Because the circuit only contains 
gates from the Clifford group, we need only simulate the 
propagation of errors that occur during the circuit [i^ 
rather than store the complete state of the quantum com- 
puter. We have written our own simulator for this pur- 
pose. We simply assume the state begins in an arbitrary 
valid codeword state and is stochastically perturbed by 
errors during the circuit. The circuit is defined to suc- 
ceed if measurement of the data qubits at the end of the 
circuit in either the X or Z basis would give the correct 
outcome. Equivalent ly, a circuit is defined to have suc- 
ceeded if an errorless error-correction cycle applied to its 
output state would produce the correct state. 

We generate data in two different ways. Where n is the 
number of levels of error correction, for n < 3, instead of 
applying an error at each location with some probabil- 
ity, we simulate the full concatenated circuit with exactly 
i errors placed randomly. For all single-qubit locations 
(state preparation, memory, measurement) the error is a 
randomly selected Pauli error and for all two-qubit lo- 
cations (CNOT, swap) the error is a randomly selected 
two-qubit Pauli error. This is repeated many times to 
generate the probability that the circuit fails given i er- 
rors, Ti. These conditional probabilities can be combined 
to give the failure rate of the circuit as a function of p, 

N 

^nr)p'{i-pf-\ (4) 

where N is the number of locations in the entire circuit, 
172, 11992, and 864496 for n = 1,2,3 respectively. The 
series is truncated after i = 6, 10, and 21 for n = 1,2,3 
respectively. These numbers are chosen so that the ap- 
proximation is valid in the region of the threshold. 

For n = 4, too much time is required to generate statis- 
tically significant results in this way. Instead, after each 
elementary physical location in the circuit we apply an 
error with probability p, where errors at all locations are 
independent. We study the failure rate of the circuit as 
we vary p between 8 x 10~^ and 2 x 10~^, again, simu- 
lating the full concatenated circuit. The time required to 
generate statistically significant results for n > 4 is pro- 
hibitive using our methods. In our simulator we use the 



SIMD-oriented Mersenne Twister pseudorandom number 
generator |42] . The message-passing method is simulated 
along side the error-correction circuit so that it operates 
in the same way as it would in a real quantum computer. 

Our results are summarized in Fig. [H The gradients 
of the lines in Fig. [H which are related to the minimum 
number of errors that cause the code to fail, are as ex- 
pected. Since the first two levels of error correction are 
unable to reliably correct errors, the failure rate for any 
less than three levels of error correction is always greater 
than p. Note that the lines for n = 3 are truncated at 
p > 10 This is because for p > 10 ~^ the approxima- 
tion that = for i > 21 breaks down as failure due to 
more than 21 errors becomes common. 

Our results suggest that the threshold is approximately 
equal to 10~^. For elementary physical failure rates be- 
low this threshold, we expect that arbitrarily accurate 
CSS operations can be performed efficiently given suffi- 
cient time and qubits. This threshold is well below the 
threshold for distillation of ancilla states that enable uni- 
versal computation [3, Ssj]. A failure rate much lower 
than 10~^ would, however, be required to achieve a suffi- 
ciently low encoded failure rate using a practical amount 
of resources. For example, we estimate that to achieve 
an encoded failure rate of 10~^^ using O(IO^) physical 
qubits per encoded qubit, a physical failure rate of ap- 
proximately 10~^ is necessary. Furthermore, if there is a 
non-zero probability of defective qubits in the linear ar- 
ray then the threshold will cease to exist. In contrast, a 
constant density of defective qubits in a two-dimensional 
array can be tolerated and will merely lower the thresh- 
old. 



V. CONCLUSIONS AND FURTHER WORK 

That the accuracy threshold for universal quantum 
computation in a system that permits only nearest- 
neighbor interactions in one dimension may be 10~^ or 
higher is somewhat surprising. In two dimensions the 
highest proven threshold for concatenated error correc- 
tion is 2 X 10~^ [lO]. There is strong evidence that the 
threshold in two dimensions can be as high as 7x 10~^ [l^ 
but this relies on techniques that are not expected to be 
useful in only one dimension. In quasi one-dimensional 
settings the highest proven threshold is 2 x 10~^ [22[. 
Since the threshold presented in this paper is based on 
numerical simulations, it would be useful to obtain a rig- 
orous bound on its value. Due to the unconventional 
error-correction method that we use, it is unclear if this 
can be done using the established level-reduction proce- 
dure of Ref. [HI. 

By adding SWAP gates where necessary, any quan- 
tum algorithm can be implemented on a linear nearest- 
neighbor array with an overhead that is, at most, poly- 
nomial in the number of qubits. This has been done 
explicitly for Shor's algorithm in Ref. [i^. Our 
result strengthens the notion that one-dimensional ar- 
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FIG. 4: Logical failure rate as a function of physical failure rate for up to four levels of error correction (lines labeled 1—4). 
No error correction (dashed line labeled 0) is shown for reference. Lines for n = 1, 2, 3 were calculated using the expansion in 
Eq. 4. To show the statistical error in the parameters there are two lines for n — 1, 2,3. The top line in each pair is 2a 
above the mean and the bottom is 2a below the mean. Data for n — A was obtained directly. The line for n — A joins the 
means of the three data points, where the error bars indicate ±2(7 statistical error. The threshold is the value of p at which 
the failure rate of the n*^ level of error correction intersects the physical failure rate in the limit n ^ oo. Assuming that the 
Unes corresponding to higher- level failure rates continue the trend observed, the threshold will be close to the crossing of lines 
for n = 3 and n = 4 at approximately 10~^. 



chitectures are viable candidates for quantura coraput- 
ing, although how to achieve defect tolerance remains an 
open question. It will be interesting to see if a higher 
threshold can be achieved by adapting the postselection 
scheme of Ref. [40| to a system that permits only nearest- 
neighbor interactions in one dimension. In this scheme 
ancilla states are protected by an error-detection code 
and postselected. We expect that the new methods of 
message passing presented in Ref. [s^ will help in im- 
proving the efficiency of such a scheme, as the weights 
outputted from error correction may be useful in moder- 
ating the amount of postselection. These ideas are the 



subjects of further work. 
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